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Abstract 


Auroral optical imagery shows marked differences between auroral features 
of the evening and morning sectors: the separation between diffuse and dis- 
crete auroras in the evening sector is not distinct in the morning sector, 
which is dominated by auroral patches and multiple banded structures aligned 
along some direction. Plasma distribution function signatures also show 
marked differences: downward electron beams and inverted-V signatures prefer 
the evening sector, while the electron spectra on the morning sector are 
similar to the diffuse aurora. We have constructed a. theory of morningside 
auroras consistent with these features. The theory is based on modulation of 
the growth rates of electron cyclotron waves by the mirror instability, which 
is in turn driven by inwar d-convected ions that have become anisotropic. This 
modulation produces alternating bands of enhanced and reduced electron precip- 
itation which approximate the observed multiple auroral bands and patches of 
the morning sector. 
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INTRODUCTION 


A major advance in the understanding of auroral arc formation processes 
in the evening sector has been achieved in the past few years by satellite 
observations and theoretical interpretations of auroral plasmas and fields 
[see, for example, Akasofu and Kan, 1981; Mozer et al. , 1980; Fennell et al., 
1981; Chiu et al., 1982; and references therein]. In particular, the electro- 
dynamic coupling between hot magnetospheric and cold ionospheric plasmas gives 
rise to a magnetic field-aligned electric potential drop which accelerates 
electrons downward to form the inverted-V and arc structures. While obviously 
a great deal of observational and theoretical details need to be examined 
further, an outline of a theoretical understanding of the relationship between 
kinetic plasma processes and the global structure of magnetospheric and iono- 
spheric plasmas and electric fields, which define the scale of the evening 
arcs, has begun to take shape IChiu et al., 1979; Lyons, 1980, 1981; Chiu and 
Cornwall, 1981; Chiu et al., 1981a, b; Schulz et al., 1981], The essential 
feature of this line of thought is schematically illustrated in Figure 1 
[Chiu et al., 1979]. The discrete-auroral oval seems to be a mapping (along 
magnetic-field lines) of a boundary layer at (or just inside) the magneto- 
spheric surface. The magnetospheric surface consists of the magnetopause and 
"neutral" sheet. The cross-magnetospheric "convection" electric field neces- 
sarily undergoes a sharp reversal (a near-discontinuity) in its meridional 
component in this boundary layer and across magnetic-field lines that map to 
the boundary layer. Note the electric field signatures on the separatrices 
between closed and open magnetic field regimes in Fig. 1. The cold colli- 
sional ionospheric plasma near the foot of the auroral flux tube is unable to 
support such a near-discontinuity in the normal (to B) component of an elec- 


1 



/ 



Fig. 1. A schematic illustration of magnetospheric structure showing 
solid closed and open field lines and their separatrices L . 
The directions of the magnetospheric electric field are shown 
by arrowheads attached to field lines. Divergence and 
convergence of electric fields occur on L . 
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trostatic field without the concomitant appearance of an electric-field compo- 
nent parallel to B. This inability is a consequence of current conservation, 
Pedersen conductivity, and ohmic resistivity in the ionosphere [Chiu, 1974 J. 
Moreover, the hoc collisionless magnetospheric plasma in the same flux tube 
supports such a parallel (to B) electric field and thus causes the E^ to 
remain substantial, even at altitudes up to about 10^ km. The signature of 
this Ej (upward in the PM sector) is plainly evident in ion and electron 
distributions observed from the S3-3 satellite [Mizera and Fennell, 1977 J. 
The indicated E^ signature at the boundary layer in the evening sector is also 
confirmed [Mozer et al., 1980 J. 

On the morningside, these conditions favorable to discrete arc formation 
by magnetosphere-ionosphere coupling are not met. Indeed, if an Eg can be 
supported in the morning sector in this admittedly simple scheme, it would be 
downward and thus accelerate electrons upward. Yet, auroral arcs due to 
electron precipitation are plainly evident in the morning sector; thus, a 
global kinetic theory of magnetosphere-ionosphere coupling remains incomplete 
unless a theory of morningside auroral arc formation consistent with the 
evening sector is formulated. 

We have constructed a theory of morningside auroras based on modulation 
of the diffuse auroral electron precipitation into bands or striations of 
characteristic perpendicular scale, which is typically of the order of several 
auroral ion gyroradii. Modulation of the diffuse auroral electron precipita- 
tion is presumed to result from modulation of the wave intensities responsible 
for scattering the auroral particles into the loss cone. This can happen, for 
example, if an ULF instability modulates the plasma density and/or magnetic- 
field strength near the magnetospheric equator and thereby modulates the 
growth rate for some VLF wave mode, either electromagnetic [Coroniti and 
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Kennel, 1970 J or electrostatic. The result could be a pulsating aurora. Such 
modulation could also occur near zero frequency and be nonconvective if the 
background plasma density is nearly uniform, as in the diamagnetic mirror 
instability le.g., Scarf et al., 1967 J. In this case, the auroral modulation 
need not be associated with micropulsations. Hasegawa [1969] has, of course, 
pointed out that the mirror instability becomes convective when the magnetic 
field is inhomogeneous, and it seems clear that a perpendicular electric field 
would also make the instability convective. However, the least convective 
excitations of the mirror instability would seem to be those for which the 
propagation vector is most nearly perpendicular to the drift velocity of the 
hot plasma, i.e. , of the plasma that is making the mirror mode unstable. The 
result should be a family of striations, roughly parallel to the discrete- 
auroral oval, but in the diffuse aurora. These striations would be passed on 
in the usual way from (a) striations in the equatorial B-field intensity and 
hot-plasma density, to (b) striations in the growth rates of electron- 
cyclotron waves (either electromagnetic or electrostatic), to (c) striations 
in the wave intensities, to (d) striations in the electron-precipitation rate 
and (perhaps more importantly) in the resonant-electron energy required for 
precipitation, and to (e) striations in the optical characteristics of the 
diffuse aurora. Such a mechanism might account for the appearance of multiple 
"arcs'' and patches in the diffuse aurora of the post-midnight sector. The 
supply of electrons for precipitation would not necessarily be depleted by 
this process, since the drift shells of hot electrons and hot protons 
intersect obliquely with each other and with the drift shells of cold plasma, 
along which the striations would tend to be aligned approximately. 

In Secton II of this report, we shall first review the plasma and optical 
observational evidences on the morningside aurora. In Section III we shall 
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discuss the properties of magnetospheric plasma distributions as they are 
acted upon by the convection electric field. We demonstrate, with interpre- 
tive support from subauroral plasma observations of the SCATHA satellite, that 
auroral ion distributions driven anisotropic by the convection electric field 
are favorable to driving the mirror instability in the morning sector, which 
is discussed in Section IV. In Section V, we offer plausible conditions under 
which a selection of the mirror instability wave number vector (k) can develop 
in the presence of shear in the convection flow. If verified, this conjecture 
would clarify the relationship between multiple auroral bands and auroral 
patches — both dominant features on the morningisde. In Section VI, we 
complete the theory by presenting the result of a phenomenological model of 
the spatial modulation of auroral electron precipitation by the mirror insta- 
bility. 

This theory of the morningside aurora, when merged with the kinetic 
theory of evening discrete arcs, may be construed as steps towards a global 
kinetic theory of auroral arc formation. 
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II. MORNINGSIDE AURORAS 


When observed on a coarse scale, the auroral oval seems to be occupied by 
a single continuous belt of auroral luminosity. However, when observed with 
finer resolution, the optical morphology of auroral arcs shows very distinct 
differences between the evening and morning sectors [see, for example, Akasofu, 
1976, 1981; Akasofu and Kan, 1981; Oguti, 1981 J. Figure 2 shows a schematic 
diagram [Akasofu and Kan, 1981 J which distills the main optical character- 
istics of the auroral morphology in various local time sectors. 

In the morning sector, Fig. 2 shows a multiple system of thin, diffuse 
and convoluted arcs in which the marked distinction between the discrete and 
the diffuse auroras is not readily apparent. Indeed, it appears that the 
diffuse aurora on the morningside is modulated into striations and patches 
[Akasofu, 1974J. A significant optical characteristic of the morningside 
aurora is that these features are frequently associated with pulsation [Oguti, 
1981; Royrvik and Davis, 1977; Johnstone, 1 978 J . Further, it has been known 
for a long time that morningside pulsating auroras have an electron spectrum 
showing a high-energy tail (> 10 keV) rather than a mono-energetic beam (< 10 
keV). This difference is illustrated in Fig. 3. Moreover, it has been 
observed [Oguti, 1976; Johnstone, 1978 J that pulsating auroral patches drift 
mostly along electric equipotentials . Thus, the pulsating auroral structure 
is defined by the warm plasma (< 10 keV), which drifts along electric equipo- 
tentials, and is not defined by the drift of the precipitating hot electrons 
(> 10 keV), which excite- the optical signature. For the hot electrons, azi- 
muthal gradient drift effects can be more dominant; thus their drift tra- 
jectories can depart drastically from convection equipotentials. The signifi- 
cant implication of this feature of morningside auroras is that (a) the forma- 
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Auroral morphology (Courtesy of Akasofu and Kan, 1981). 
Features of Interest to this report are multiple striations 
and patches of the morningside. High-latitude arc systems are 
to be distinguished from these. 
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cion mechanism of the striations and patches of auroral luminosity must be 
sought in irregularities of the warm (< 10 keV) plasma, and (b) the irregu- 
larities of the warm (< 10 keV) plasma must be able to modulate the precipita- 
tion of hotter electrons which precipitate and produce the auroral luminosi- 
ty. We shall see in the following that our theory is consistent with these two 
observational implications. 

Although morphological features cannot be discerned from satellite plasma 
observations, satellite plasma data on the morningside are a powerful indica- 
tor of the kinetic state of the plasma forming the aurora. Unfortunately, due 
to the predominant research interest in evening inverted-V structures, very 
little auroral plasma data interpretation for the morningside has taken 
place. Figure 3a shows a typical evening auroral electron spectrum measured 
by the S3-3 satellite. In contrast to the electron beam features character- 
izing the evening discrete arcs of Fig. 3a, the morningside electron spectrum, 
shown in Fig. 3b, consists of a hardening at ring current energies and no 
mono-energetic bean. This is entirely consistent with modulation of hot 
electron precipitation as the cause of multiple morningside arcs. Although no 
formal statistical measure of the local time preference of these two types of 
electron spectra has been obtained, the predominant preference of inverted-V 
precipitation events for the evening sector is well-established. For example, 
this tendency is plainly evident in the AE-D data [Hoffman and Lin, 1981 J. 
Further, a statistical survey of upgoing ion beams, another element of plasma 
signature of inverted-V structures, has been made [Gorney et al., 1981] and 
the results show definite preference for ion beams to occur in the evening 
sector while few ion beams are found in the morning sector. A particularly 
succinct summary of the differences between auroral plasma characteristics 
at ~ 1 Rg altitude is shown in Fig. 4. In this symmetrical pass of the S3-3 
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S3 3 PARTICLE WAVE SPECTROGRAM JUNE Ifl, 1977 
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Fig. 4. Energy-time spectra of a dawn-dusk pass of the S3-3 

satellite. The top panel is a wave frequency vs. time plot of 
the pass. 
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satellite across the auroral oval, it is seen that the simultaneous occurrence 
of an electron inverted-V and an ion beam indicates parallel potential drops 
in the evening sector both above and below the satellite, while, on the mor- 
ningside, the absence of both features in the plasma observations is evidence 
that the morningside aurora is not associated with electric potential struc- 
tures, either above or below the satellite. (For a detailed interpretation of 
features such as in Fig. 4, see Chiu et al., 1982.) 

At whatever local time, the plasma distributions at the equatorial 
regions of auroral field lines are crucial to any aurora formation theory. 
Unfortunately, direct ' measurements of equatorial auroral plasmas are not 
available on a continuous basis. However, invoking the concept of plasma 
convection in the magnetosphere, features of equatorial auroral plasmas can be 
reconstructed from equatorial subauroral and near-auroral observations such as 
those of the SCATHA satellite [Fennell et al., 1981 J. Since these data are 
crucial to the construction of our theory, we need to devote special attention 
to them; therefore, we defer discussion to the next section where convection 
of plasma distributions in the magnetosphere is treated. 

Meanwhile, we note the following observational features which distinguish 
morningside from eveningside auroras: 

1. Multiple thin arc systems and patches dominate in the morning, and 
single discrete arc in the evening. 

2. Inverted-V electron precipitation structures and ion beams are 
primarily evening features, indicating absence of parallel electric accelera- 
tion in the morning. 

3. Morningside electron spectra, similar to the typical diffuse aurora, 
generally show power-law type hardening at energies higher than a few keV. 
These spectra are distinctly different in that the mono-energetic downward 
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electron beam is missing. 

4. Morningside optical features are seen to drift with the warm (< 10 
keV) plasma, not with the energetic electrons which produce the optical fea- 
tures themselves . 

5. Pulsating patches, a predominantly morningside phenomenon, exhibit 
characteristics consistent with items 3 and 4 above. 

6. Quasi-periodic recurrent auroral patches exist in the dawn sector. 
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III. FORMATION OF PLASMA DISTRIBUTIONS 


To construct a global theory of auroral arc formation, one must incorpo- 
rate the plasma and field features dictated by global boundary conditions, 
which in turn provide information about how and from where the hot auroral 
plasmas originate. Virtually all current ideas on the supply and maintenance 
of hot auroral plasmas invoke inward convection of plasmasheet particles. In 
the simplest form, the inward convection hypothesis assumes that a single 
uniform population of isotropic ions (T^ ~ 2 - 3 keV) and electrons (T g ~ — - 
1 keV), found well inside the plasmasheet (Frank et al. , 1981], is convected 
adiabatically inward past the auroral zone and perhaps deep into the plasma- 
sphere. This simple theoretical concept is appealing in many respects, but it 
requires testing by comparison with plasma observations in (or near) the 
equatorial regions of auroral and/or sub-auroral field lines. We shall see in 
this section that plasma observations of the SCATHA satellite imply essential 
modifications of this simple concept. These modifications of the inward 
convection hypothesis of hot plasmas, taking into account hot ions of iono- 
spheric origin, are a key element in constructing our theory of morningside 
multiple auroral bands, although the most important element of time-dependence 
[Chiu and Kishi, 1983] will be dealt with elsewhere. 

When the SCATHA satellite was recognized to be near the equatorial 
regions of auroral (and/or near auroral) field lines, the plasma instruments 
onboard measured typically two types of ion distribution functions and a 
typical electron distribution function with anisotropy perpendicular to the 
magnetic field (3). The two types of ion distribution functions are shown in 
Figs. 5a and 5b, and the typical electron distribution function Is shown in 
Fig. 6. In both ion and electron distributions, it is seen that the more 
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energetic populations are anisotropic in the v ^-direction, as one would expect 
if these particles were convected in from the plasmasheet. For the ions, a 
low-energy component, anisotropic in the v^-direction, is almost always 
present but in varying strength (see below). The presence of 0 + in this 
field-aligned component indicates strongly that it originates from the dis- 
crete aurora [Fennell et al., 1981; Kaye et al,, 1981 J. The field-alignment 
and composition strongly suggest that these ions are the low-energy remnants 
of ion beams accelerated upward by the electric potential structure associated 
with the auroral discrete arc. Since the discrete arc and inverted-V struc- 
tures are primarily evening and midnight phenomena, the relative dominance of 
these low-energy field-aligned ions should also reveal a preference for the 
evening and midnight sectors. Fennell et al., [ 1 98 1 J distinguished the field- 
aligned dominant distribution of Fig. 5a from the field-perpendicular dominant 
distribution of Fig. 5b by the value of the energy corresponding to the separ- 
atrix between the field-perpendicular and field-aligned distribution function 
contours (the transition energy). Thus, the transition energy of Fig. 5a is 
greater than that of Fig. 5b. Fennell et al., [1981] showed a plot of this 
transition energy as a function of magnetic local time for almost three days 
of SCATHA data [Fig. 7J. Although real-time geomagnetic activity obviously 
plays a role in the increase of the field-aligned population, we note that 
increases generally take place in the 18-03 MLT sector where the discrete arcs 
are generally located. Further, as one approaches the morning sector (after 
04 MLT), the field-aligned population shows a fading trend. Thus, at the 
morning sector, on near-auroral field lines, the dominant energetic ion popu- 
lation is that of plasmasheet origin, as indicated by the perpendicular pitch 
angle anisotropy of Fig. 5a. Although the demarcations between "evening-like" 
and "morning-like” trends in Fig. 7 are not sharp and vary from day to day, we 
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Fig. 7. Local-time dependence of the transition energy between field- 
aligned and field-perpendicular populations such as shown on 
Figs. 5a and 5b. 
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stress that these trends are consistent with the preferred locations of ion 
beams and inverted-V's, and with the inward convection hypothesis. In this 
report we shall assume the above interpretation and examine its consequences. 

We shall argue in the next section that the field-perpendicular popula- 
tions of ions favor generation of mirror instability at or near zero-frequency 
in the morningside auroral zone, whereas the presence of the lower-energy (< 
10 keV) field-aligned population in the late evening sector suppresses mirror 
instability in the auroral region, thus accounting for an important difference 
between evening and morning auroras. 

From the above description and qualitative interpretation of the distri- 
bution functions and transition energies of Figs. 5-7, it is rather obvious 
that the appearance of hot plasmas in auroral and sub-auroral equatorial 
regions can no longer be considered in terms of time-independent convection of 
a single component of plasma from the plasmasheet. In short, the causal 
implications of injection in real time, the local -time dependence of the 
injection source, and, last but not least, the existence of an ionospheric 
source from the discrete aurora are essential and unavoidable complications 
that are part and parcel of the inward convection hypothesis for the formation 
of hot auroral plasmas [Chiu and Kishi, 1983J. Thus, in order to verify 
theoretically the interpretation of data discussed above, we need to construct 
nothing less than a causal theory of plasma injection and convection in the 
equatorial magnetosphere, taking into account spatial and temporal factors of 
both the sources and the convection of plasma distributions. A Green’s func- 
tion solution of such a new time-dependent theory has been formulated and will 
be dealt with elsewhere [Chiu and Kishi, 1983]. For our purposes here of 
interpreting Figs. 5a and 5b for the construction of auroral ion and electron 
distributions, it is sufficient to abstract the local-time (spatial) trend of 
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Fig. 7 and Co focus upon Che spaCial source discribucion underlying Fig. 7, 
racher Chan co focus upon Che time-dependence of Che injecCion process which 
complicates Fig. 7 in Che midnighC hours. Under Chis simplif icaCion, morning- 
side and eveningside ion discribucions can be considered in cerms of a time- 
independenc plasma convecCion model, which is Che presenc "sCaCe of Che arc" 
le.g. Wolf, 1970; Chen, 1970; Kivelson eC al., 1979; Kivelson and SouChwood, 
1975; Cowley and Ashour-Abdalla, 1975; Kaye eC al., 1979 J • 

In such a sCaCic model, plasma is driven by Che conveccion-corocacion 
eleccroscacic poCenCial [Volland, 1975 j 

4 (L, - $ c L 0 * 1 sin* + $ r /L (1) 

where n=l is an assumed parameCer, 4> c is a measure of Che cross-Cail poCenCial 
drop, and 4^ ■ -91 kV is a measure of Che dipolar corocacion poCenCial. Figure 
8 shows an equipoCenCial concour ploc of (1), wich 4^ ■ -0.43 kV (roughly 
appropriaCe for storm time). For the plasmasheet source distribution, we 
assume a uniform nightside (dawn to dusk) Maxwellian distribution of ions and 
electrons at L“13. In Che evening sector, we assume field-aligned ion beams 
are to be injected at L=10 by boch the northern and southern discrete arcs, 
since distinct counter-streaming ion beams are often seen even at SCAThA 
latitudes. The ion beams are assumed to have drift energies of < 10 keV and 
Maxwellian temperatures of similar magnitude. Plasma distributions from these 
two sources are assumed to- be adiabatically convected to the location of 
SCATHA at L ~ 7 at various local times. The source distribution function is 
given below: 


f - N m [n^/2* T m J 3/2 • exp[-m^ v 2 /2T M J • g(L, 4) 
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Equipotential contours (in kV) of a simple convection electric 
field model according to Eq. (1) (Volland, 1975). 
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where the subscripts M and A denote plasmasheet and discrete auroral sources, 
respectively. The spatial functions g and h are as described above: 


g - <S(L - L m ) 9 (<fr -v/2) 9 (3w/2 - <P) 


(3) 


h - fi(L - L ) 9($ - v/2) 6 (it - <*>) 


(4) 


where $ is measured clockwise from noon MLT ($ ■ 0). 

The adiabatic mapping of a particle of mass m and charge qe from the 
phase space point (Vq, Tq) = (vq, a^, Lg, to the phase space point 

(v, r) = (v, a, L, <P) requires conservation of the magnetic moment 

v L 3 '' 2 sin a - v Q Lq 3 ^ 2 sin a Q (5) 

and conservation of energy 



qe4(L 0 , <J»q) 



+ qe $(L, <p) = W 


( 6 ) 


In addition, the conservation of the bounce invariant is assumed to be given 
by the approximate relation [Southwood and Kivelson, 1975] 


sin a , . , , T /t n 0.45 ,, . ,,-l 

iI — = [sin a Q + (L/L 0 ) (1 - sin a Q )J 


(7) 
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Figures 9a and 9b show an example of our source model distribution napped 
to two observation points: (a) L * 7 , $ * tf/2 and (b) L ■ 7, ■ 3r/2. We 

have assumed T^ = 1 keV and T A ■ 1/2 m A v^ » 0.5 keV. These model 
parameters, together with the relative values of and N A , can be chosen to 
fit the important population components of Fig. 5a. The upper energy cutoffs 
in Figs. 5a and 9a are the effects of Alfven boundary formation at the L-value 
of observation. The field-alignment of the lower-energy population in Fig. 9a 
is due to the merging of two counter-streaming beams of the model, (4). 
Although Fig. 5a does not show counter-streaming beams at low energies, the 
counter-streaming ion beam phenomenon is frequently observed in SCATHA data 
[Richardson et al., 1981 j • Since the signature of counter-streaming ion beams 
is not always present but field-aligned anisotropy is always present in the 
low-energy component, the auroral source in (4), consistent with Figs. 5a and 
5b, can also be modeled as a bi-Maxwellian with field-aligned anisotropy. 
Indeed, it seems reasonable that the two-stream instability could cause two 
counter-streaming beams to merge into a field-aligned population [Schulz and 
Koons, 1972 J . This simpler bi-Maxwellian source model, with proper convection- 
induced cutoffs, is used in illustrating mirror instability calculations, 
although both counter-streaming and bi-Maxwellian source models can be used 
with similar effect. 

To conclude this section, we note that the systematic local-time related 
features of SCATKA plasma populations observed at L ~ 7 can be understood in 
terms of relatively simple sources at higher L-values together with the simple 
concept of inward adiabatic convection. This exercise gives us confidence in 
the model for L-values between L = 7 and the sources. 

We shall demonstrate in the next section that morningside field- 
perpendicular ion distributions such as Fig. 9b can drive mirror instability 
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MLT = G hr 


MLT = 10 hr 




A. MORNINGSIDE n. EVENINGSIDE 


Fig. 9 . Simple Maxwellian ion distributions mapped from L = 13 to (a) 

L = 7 and dawn, and (b) L = 7 and dusk. Field-aligned ion 
beams are assumed to be injected at T, = 10 in the evening 
hours . 




at L 3 7-9, whereas the field-aligned component in the evening, such as in 
Fig. 9a, can stabilize ' the effect of the field-perpendicular component so as 
to inhibit mirror instability in the evening hours, if the auroral source can 
provide sufficient field-aligned ions of sufficiently high energy. 
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IV. MIRROR INSTABILITY AND PLASMA CONVECTION 


It is known for a long time that a high-3 (3 ~ 1) plasma with pressure 
anisotropy (3^ > 3^) can be unstable to a longitudinal wave mode which is an 
offshoot of the MHD waves. This instability was named the "mirror" instabili- 
ty because a loss-cone distribution in a mirroring field would have 3^ > 3^ . 
In space plasma physics, the mirror instability has frequently been invoked to 
account for compressional mode micropulsations [Hasegawa, 1969; Lanzerotti et 
al., 1969J with periods in the > 100 sec range. Hasegawa [1969] showed that 
in a plasma density gradient the mirror mode merges with drift waves such that 
the real part of the frequency is the drift wave frequency. At lower 
L-values , such as near the plasmapause region, such micropulsations, with non- 
zero frequency and interpreted in terms of the mirror instability, may indeed 
be very Important. Since we are interested in the mirror instability driven 
by convected plasmasheet ions In the region poleward of the diffuse auroral 
boundary (mapped to the equatorial plane) in the morning sector where there 
are no marked cold plasma gradients, we shall not be involved with the micro- 
pulsation aspects of the mirror instability, although particle flux modulation 
is part and parcel of micropulsations. We emphasize that zero -frequency 
mirror instability need not be correlated to occurrence of micropulsations. 
In the absence of the real part of the frequency, the mirror Instability 
becomes striations of the plasma and furrows of the magnetic field in space. 
In the next sections, we shall demonstrate that these striations of the con- 
vecting plasma and magnetic field can affect the rate of diffuse auroral 
electron precipitation. 

To determine whether convecting ions can drive the mirror instability, 
one must show that the distribution functions modeled in the last section 
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satisfy a criterion [Hasegawa, 1975, pp. 94-99J under which the growth rate 
can be positive. In the case of bi-Maxwellian plasma, the instability cri- 
terion is 

1 + 3 ii (1 " T li /T Ji } + S le (1 " T le /T Ie ) < 0 (8) 

2 

where 0^ = 8nnTj/BQ for both ions (i) and electrons (e). It can be seen in 
(8) that > Tj and 0^ ~ 1 are required to satisfy the instability cri- 

terion. It is our contention that both conditions are satisfied by plasma- 
sheet ions and electrons injected onto the auroral region because inward- 
convecting particles gain perpendicular energy preferentially, and substorm 
plasma at high L-valued equatorial regions are expected to have higher 0. In 
order to demonstrate instability in our model of convected plasma distribu- 
tions, (8) cannot be used because the distributions shown In Figs. 9a and 9b 
are obviously not bi-Maxwellians , even though some sources can be bi- 
Maxwellian. 

In the notation of Hasegawa [1969J, in which e is the i-j component of 
the plasma dielectric tensor, the relevant dispersion relation to mirror 
instability in the low-frequency (w/fl^ < 1) and long-wavelength < 1) 

regime is 


k 2 - ^ (e + e 2 /e ) 

c 2 yy yz zz 


(9) 


In the above, k is the wave number, w is the wave frequency, is the cyclo- 
tron frequency of species j = i,e, and Vj is the thermal velocity. The third 

2 2 2 2 

term of (9) in the form (w /c k )(e /e ) is of the order of the 

yz zz 

electron 0 g (< 1); therefore, the dispersion relation is assumed to be k“ 
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- c/e /c = 0, where in terns of a general distribution function f. and 

yy J 


plasma frequency cu 


a =i+z (u 2 ./n 2 ) + 

yy jai>e pi j 


£ 21 i dv i l v i dv i ( k i > v ? - i v ! J £ j - k i v i f j 

j-i,e w J -* 0 j j 


.2 1 2 (5 V 3 V, 

- k ! • 4 V 1 v,- 4,yic, ) 
B H 


( 10 ) 


From (9) and (10), it is a straightforward matter to show that the instability 
criterion is 


1+^f 2 

B q J-i,« 


(3f . /3v . ) 


2r 


I dv, / dv i (i « vj £ + i. v d P - j - a/ M < 0 an 
-« o j j o n 


where P indicates the principal value integral. In the low-frequency limit 
(indeed, ui ~ 0 in our case) the resonance denominator in (11) can be expanded, 
and (8) is obtained if fj is bi-Maxwellian. 

Since the empirical distributions (Figs. 5a and 5b) and the model distri- 
butions (Figs. 9a and 9b) are not bi-Maxwellian, a realistic evaluation of the 
instability criterion (11) requires numerical techniques. We have applied the 
convected distributions of source models discussed in the previous section to 
the instability criterion (11) to show that, if the field-aligned auroral ion 
beam source is sufficiently strong (in number and in energy) , the mirror 
instability driven by the field-perpendicular plasmasheet component can be 
stabilized in the evening sector, thus confining the effects of the mirror 
instability to the morning sector during strong injection events such as in 
Fig. 7. To illustrate this point, the source model (4) is not very convenient 
because the field-aligned anisotropy of the higher-energy portion of the 
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convected distribution from the auroral source depends on T^ and v A in a 
complicated manner. In contrast, the source model with a bi-Maxwellian 
auroral source (with temperatures and T^) is particularly convenient 
because the field-aligned anisotropy of the convected auroral component (Fig. 
9b) is controlled by the source field-aligned anisotropy ^^/T^. Figure 10 
illustrates the dependence of the indicator of mirror instability, left hand 
side of the inequality (11), as function of at the source for a given 
normalization density ratio between the plasmasheet and auroral components of 
Fig. 9b. Thus, if the discrete aurora is sufficiently strong and the ion beam 
energy is sufficiently high (i.e., the source ions are strongly field- 
aligned), the mirror instability can be stabilized by the evening discrete 
arc. On the other hand, the absence of the discrete auroral source on the 
morningside enhances the likelihood of mirror instability associated with 
injection events of Fig. 7. For the calculations of Fig. 10, the magneto- 
spheric source is located at " 11 instead of % = 13. For = 13, insta- 
bility can occur at even higher values of T j/T^ than that of Fig. 10. 

This new concept of "feedback" by the discrete aurora upon the effects 
driven by magnetospheric convection is somewhat interesting in that it breaks 
away from the perception that the aurora is a passive sink driven by solar 
wind-magnetosphere interactions. 
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MIRROR INSTABILITY INDICATOR AT L 



Fig. 10. Dependence of Che mirror instability indicator at L, = 9 [left- 
hand side of the inequality (11)] on the ion anisotropy T f /T^ 
of the discrete aurora. L 
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V. SELECTION OF £ VECTOR 

In order to construct a theory of the multiple bands and patches of the 
morning sector, it is not sufficient to merely establish instability. A 
plausible explanation of the scales of the bandedness is necessary. This 
means examining mechanisms by which a single typical value and direction of 
the instability wave vector k are selected. 

For a given k the mirror instability manifests itself as a series of 
growing rarefactions and condensations of magnetic field and plasma occurring 
in parallel planes. However, the disturbance thus represented by a given k 
does not propagate as if in an elastic medium. Rather, it moves only in 
association with the drift motion of the underlying plasma. This is the 
significance of having Re u> *» 0 in (8) - (II), where the mirror instability is 
analyzed in the drift frame of a uniform plasma. 

It follows that in a uniform plasma there is no preferred direction for 
k. In that case the mirror instability should take the form of a jumbled 
superposition of k-vector directions and thus resemble (at least in the plane 
perpendicular to B) the condensation of droplets from a vapor. This morpholo- 
gy may be relevant to the patchiness seen in the equatorward portion of the 
diffuse aurora in the 04-13 sector of local time (see Figure 2) but cannot 
describe the azimuthally elongated bands that seem to evolve from the entire 
diffuse aurora soon after midnight and to persist in the poleward portion of 
the aurora throughout the AM sector. The presence of such a banded structure 
suggests the preferential selection of k vectors perpendicular to the plasma 
drift velocity. 
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A likely selection mechanism involves the presence of velocity shear in 
the magnetospheric plasma motion. The effect of velocity shear on the mirror 
instability is to refract the wave-normal vector k, with respect to the 
plasma-velocity vector u, so as eventually to align the wavefronts with the 
streamlines. The refraction process is analogous to the rotation of ranks 
that would occur down a broad avenue if a phalanx of soldiers were to march in 
the x-direction with speed u and perpendicular velocity shear du/dy. Let k 
be the rank normal vector. The separation of ranks along the x axis remains 
constant at since the velocity of march is independent of x. The slope 

of a line normal to the ranks is ± ky/kjj. The ambiguity of sign stems from 
the possibility of constructing either a forward (+) or a backward (-) rank 
normal. It follows that for rectilinear motion the rank-normal direction 0 
evolves according to 

tan 9(t) = k /k * tan 9(0) + (du/dy)t (12) 

y x 

and approaches + ir/-2 as t + + “, as. one would expect. 

The ranks of soldiers in this allegory represent plasma condensations or 
rarefactions associated with the diamagnetic mirror instability, and the 
marching velocity u represents the plasma-drift velocity across the magnetic 
field B. The rank-normal vector k represents the wave-normal vector, and the 

A A 

rank-normal angle 9 represents the wave-normal angle cos -1 (k • u). The 
foregoing results already imply that 9 approaches + it / 2 as t + + •, i.e., that 
the plasma condensations and rarefactions become aligned with the drift- 
velocity streamlines. However, in the plasma context it is equally important 
to note that 
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9 9 1/2 9 9 1/2 

k.(t) = (k“ + k‘) = { [k (0)T + [k v (t)r} 

-l x y x y 

= {l + [tan 6(0) + (du/dy)t]^} |k x (Q) | (13) 

where k^(t) is the magnitude of the vector k * B. For * 0 one finds that 

k^(t) will increase without limit as t -*• + “. Now since the maximum growth 

occurs for k^P^ ~ 1, where P i is the mean ionic gyroradius [Hasegawa, 1969 J , 

modes with * 0 are not favored for growth. Thus, the presence of shear in 

the plasma drift-velocity field favors the excitation of mirror-mode stria- 

tions for which k„ is initially of order p. ^ and for which k, is initially 

y i A 

zero. The striations, in other words, are preferentially aligned with stream- 
lines of the velocity field. 

The magnetospheric situation is slightly complicated by the fact that 
"hot magnetospheric plasma" does not have a unique drift velocity at any point 
in space, since the gradient-curvature component of that drift depends on the 
energy and pitch angle of the individual particle. A unique drift velocity 
could be defined, of course, by averaging over the energy spectrum and pitch- 
angle distribution at each point in space. However, it should be sufficient 

-V 

for the purpose at hand to observe that shear is present in the E * B drift 
associated with magnetospheric convection, and that gradient-curvature drifts 

■f «*■ 

are typically much smaller than E * B drifts for particles representative of 
the inner plasmasheet. For example, the ratio of gradient-curvature drift to 
convection velocity at L ~ 8 is roughly the ratio of particle energy to 10 
MeV. 

The magnetospheric situation is further complicated by the fact that 
E x 3 drifts are not rectilinear, thus a tensor quantity is needed to express 
the presence of shear in the drift-velocity field. One must analyze the shear 
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tensor [e.g., Morse and Ingard, 1968, pp. 272-273J specified in rectangular 
coordinates by the expression 

S t * (1/2) t(3u 1 /3xj) + Ouj/te^J - (6^/3) 7 • u (14) 

in which x^ = x, X£ = y, and x^ ■ z. For the field model illustrated in Fig. 
7, both E and B are derivable from scalar potentials. It follows that 

7 • u - 7 • [(c/B 2 )(E x f)J 

= (c/B 2 ) [B • (7 x E) - E • (7 x B)J 

- 2 (c/B 3 )(E x b) • 7B 

=> - (2/B) n • 7B (15) 

after application of well-known vector identities [e.g., Panofsky and 
Phillips, 1962, p. 471 J. The electric field E(= -70) is derivable from (1) if 
the dipolar relationship L = (r/a) csc^ 0 is inserted there. The coordinates 
r and 9 denote geocentric distance and colatitude, respectively. The magnetic 
field is given by 


B - (p/a 2 ) 7 l(a/r) 2 cos 9J 


(16) 


where a and p represent the radius and magnetic moment, respectively, of the 
earth. Evaluation of u = (c/3 ) E x B in the equatorial plane leads to the 
simple expressions 
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u = - (c/u) [(4 /a 2 )(x 2 T y 2 )(x 2 + 2y 2 ) - 4 ayj 
x c r 


(17a) 


and 


u » (c/y) [(4 /a 2 )(x 2 + y 2 ) xy - 4 ax] (17b) 

Y c r 

for the plasma-velocity vector u. The nonvanishing components of the shear 
tensor in this case are 

S u =* - 2(c/y)(4 c /a 2 )(x 2 + 2y 2 ) x (18a) 

S 2 2 = ( c/y)($ c /a 2 )(3x 2 + 5y 2 ) x (18b) 

S33 = - (c/y)(4 c /a 2 )(x 2 + y 2 ) x (18c) 

and 

S 12 “ S 21 = - (c/2y)(4» c /a 2 )(3x 2 + 7y 2 ) y (18d) 

where x = r cos <j> and y = r sin <j> in the equatorial plane. The calculation of 

S33 in (18c) requires slightly devious methods, since it is not true that 

3u z /3z 3 0 in the equatorial plane. However, it follows from (17) and (15) 

that (3u /3x) + (3u /3y) = (1/2) V • u, and therefore that 3u /3z = (1/2) 

x y z 

7 • uj this implies that S33 * (1/6) V • u, and inspection of (18) reveals 
that the shear matrix is indeed traceless, as is required by (14). Note that 
(18) does not contain 4 , the parameter associated with rigid rotation of the 
magnetosphere . 

Of special interest in the present context is the quantity 
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s = 


(19) 


2 

ij 


“i s ij e j 


2 

ij 


S. . 
ij 


u . 
J 


which represents the velocity shear measured with respect to a displacement 
normal to the velocity itself in the plane perpendicular to B. The symbols 

a A 

and u^. in (19) denote Cartesian components of the unit vectors parallel to 
E and u, respectively. Evaluation of (19) in the equatorial plane yields 


s ■ - [ r^ (13x“ + 28y^)y - 2(4^ a 3 /4 c ) r 2 (5x 2 + 14y 2 ) 

+ 7 <# r a 3 /$ c ) 2 yj r 2 (c<> c /2a 2 y) 

+ [r^ (x 2 + 4y 2 ) - 4(4 a 3 /4 ) r 2 y + (4 a 3 /4 ) 2 J (20) 

r c r c 

where r^ = x^ + y^. Both numerator and denominator of (20) vanish at y = 
3 1/3 

(4 a /4 ) on the dawn-dusk meridian (x = 0), but the expression for s 
r c 

remains continuous and approaches the value -(7c/4y)4 r a there. This is the 
point at which equipotentials seem to cross in Fig. 7, and at which the plasma 
drift velocity therefore vanishes. 

Contours of constant s in the xy plane are a measure of the direction 
along which striations will be approximately aligned. These are shown in Fig. 
11, which is to be considered together with the equipotential contours from 
Fig. 8. The amount of shear thus seems to increase with r in the region of 
open equipotentials, and therefore in the region corresponding to the auroral 
zone. This means that the selection mechanism for k becomes increasingly 
effective with increasing L, in accordance with the observational pattern (J. 
C. Siren, personal communication, 1982) illustrated in Figure 2. In other 
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Fig. 11. Contours of constant shear s [defined in Eq. (19)] of the 
convection flow corresponding to the electric potential of 
Fig. 8. 
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words, the morningside portion of the diffuse aurora evolves into azimuthally 
elongated bands (suggesting a unique k) at the higher L values but into a 
patchy morphology (suggesting a jumbled superposition of k vectors) at the 
lower L values. Since the flow model described by (12) yields s = 
( l/2)(du/dx), which is roughly appropriate on the local scale (as distinguished 
from the global scale of (20)), we presume that the "wave-normal” angle 0(t) 
and "wave vector" k(t) for the mirror instability will evolve roughly in 
accordance with the equations 

tan 9(t) = k • E/k • u ■ tan 9^ + 2 / s(t') dt* (21a) 


and 


t rp 1 /2 ^ A 

k ± (t) = {l + [tan 8 q + 2 / s(t') dt'P} |k • u| Q (21b) 

as one follows the E * B drift of a plasma element in more realistic flow 
geometries, such as that which leads to (13). The subscript 0 in (21) denotes 
evaluation at time t * 0. 

We have so far neglected the effect of velocity shear on the plasma 
stability of the flow-velocity field itself. However, it is well known that 
velocity shear can lead to a Kelvin-Helmholtz instability whereby irregulari- 
ties and eddies develop in the velocity field. Familiar examples of this from 
everyday experience include wind-driven surface waves on a body of water and 
flags waving in a steady breeze. The effect of Kelvin-Helmholtz instability 
on auroral forms was analyzed by Hallinan and Davis [ 1970 j . Figure 2 contains 
clear evidence of this effect on the azimuthally elongated forms that evolve 
out of the diffuse aurora, as well as on the discrete arc. 
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One might question our use of a dipolar magnetic-field model in the outer 
magnetosphere, to which the auroral oval maps, for the purpose of calculating 
velocity shear. Indeed, it is not very difficult to calculate as a func- 
tion of r, 9, and $ for a model magnetosphere like that shown in Fig. 1. 
Since our essential points concerning the selection of k vector and its con- 
nection to morningside auroral morphology can be made with the simpler dipolar 
calculation, we shall not be involved with such a complex case. 
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VI. MODULATION OF CYCLOTRON INSTABILITY 


The above results suggest that the mirror instability creates spatial 
irregularities in plasma density and magnetic-field strength in the equatorial 
portion of the outer magnetosphere, and that such irregularities are typically 
aligned with the plasma-convection pattern. We presume that such irregulari- 
ties modulate the wave-particle interaction that is responsible for the dif- 
fuse auroral precipitation of energetic electrons. The precipitating 
electrons have perhaps tens of keV of energy and interact most probably with 
electrostatic electron-cyclotron waves, since the diffuse aurora occurs well 
outside the plasmasphere. The particles responsible for the mirror 
instability are mainly ions of the same energy range. In other words, the 
mirror instability modulates the electron-cyclotron instability, and 
altogether different particle populations are involved in the two 
instabilities . 

We analyze the effects of this modulation on the intensity of electro- 
static-cyclotron waves and on energetic-electron precipitation by means of the 
model equations introduced by Schulz [ 1974a. J to simulate the Kennel-Pets chek 
1 1 966 J process. These equations are of the form 


dl 

dt 


II 

I + Xx 


+ S 


(22a) 


and 


1 d\ 2 Y* (I/I*) ,i+ j in R , W 
I dt = ~ 1~T + l v g • B l IT + X 


(22b) 


and they are designed to simulate the sighificant physical phenomena without 
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imposing the heavy calculational burden that an analysis from first principles 
would entail. The symbols A and I in (22) correspond to wave and particle 
intensities, while the symbols W and S represent their respective source 
strengths. For our considerations here, we assume that all electrostatic 
cyclotron waves are generated by the energetic electrons of interest, so W » 
0. The particle lifetime approaches the strong-diffusion lifetime T (~ 100 
sec) in the limit At » 1 and approaches 1/A in the weak-diffusion limit 
At « 1 . The growth rate for wave amplitude is represented as (1 

+ At) * y (I/I*). The factor (1 + At) * implicitly describes the dependence 
of the growth rate on pitch-angle anisotropy, which vanishes in the limit of 
strong diffusion (At » 1). The normalizing factor I* for particle intensity 

is analogous to the intensity limit described by Kennel and Petschek [ 1966 J in 

* 


the context of unstable whistler-mode waves. The factor Y expresses the 

proportionality of growth rate to particle intensity at fixed anisotropy. The 

term |v * B/La| In R in (22b) describes the fact that waves have only a 
§ * 

certain probability R of survival after propagating a distance ~ La along the 

field line at group velocity v , where L is the magnetic shell parameter and a 

S 

is the radius of the earth. Thus, the parameter R serves as an effective 
internal-reflection coefficient for wave energy. 


The time derivatives that appear in (22) must be interpreted as total 
(Lagrangian) time derivatives in the reference frames of the particles and 
waves, respectively. The difference between reference frames must be recon- 
ciled before (22a) and (22b) can be solved self-consistently. Thus, we must 
write 


dl/dt - (3I/3t) + u • V ± I + r£2 3 <p • V ± I 


(23a) 
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and 


dX/dt » (3X/3t) + u • Vi + v • V.X (23b) 

1 S 1 

where u is Che convection velocity (c/B^) B x b and rft^ is the bounce- 
averaged gradient-curvature drift velocity (mapped to the equator) for the 
average energetic electron. The spatial derivatives in (23) involve 
rather than V because (22) represents (as usual) a bounce-averaged treatment 
of the wave-particle interaction. 

We seek a steady-state solution of (22), a solution for which the impor- 

* * 

tant spatial variations of X, I, and Y /I occur in the direction of the 
ic vector of the mirror instability, i.e., in the direction of the electric 
field E. The tendency of strong shear in the convection velocity u = (c/B^) 
E x b to select this orientation for k has been demonstrated above. We 

reduce the steady-state problem from two spatial dimensions to one by neglect- 

* * •+■ ■+• _ 
ing any variation of X, I, or Y /I in the direction of u, i.e., u • *» 

0. Thus, we sacrifice a description of global structure in order to concen- 
trate on spatial modulations produced locally by the mirror instability. 

Thus, we can define a local Cartesian coordinate system (<, y) such that the 
* -*■ -*■ 

unit vectors < and y are parallel to k and u respectively. In this system, 

A g a 

V = ic - 3 — and the furrows of the mirror instability lie parallel to y, as 

X oK 

illustrated in Fig. 12. Note the important difference between the gradient 

A A A 

drift direction <p and the local vectors (k, y) oriented by the shear action in 
the convection flow. 

A further notational simplification of (22) can be introduced by 
normalizing to the weakly diffusive (Xt « 1), unmodulated Kennel-Petschek 

* ic 

limit without residual wave source (W => 0), i.e., I = Iq, Y = Yq and dX/dt = 
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Fig. 12. Coordinate system used in considering electron flux 

modulation. $ is the azimuthal direction, k is the wave 
vector direction, and u is the convection flow direction. 
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0. In this limit (22b) yields 




A 

B| 


in R 
L a 


(24) 


which can be reintroduced in (22b) to eliminate the constants involved. 
Further simplification can be' achieved if I and Y are to be measured in units 

it it 


of I Q and Yq. 


Under these assumptions, the model equations become, in the frame of the 
flux tube drifting with velocity u. 


- - 3 (i/v Xt (I/ V * 

T r a 3 (♦ • K) - - -p r rr- + Cxs/I 0 ) (25a) 


_1 

\ 



3 <- Xt > , 2 y*t(I/I*) _ - * 

; 3k 1 + X t ZT 0 


(25b) 


where the source S and strong-diffusion lifetime t are regarded as con- 
stants. The periodic spatial variability of interest originates in the 

it it 

ratio (y /I ) in (25b), which we model as a plane wave modulation of the 

* * 

unmodulated Kennel-Petschek limit (Yq/Iq) 


Y*/I* ■ (Yq/Iq) exp (Q sin k. ic) 


(26) 


We expect the parameter Q in (26) to be of order unity. This estimate is 
based on the expectation that the mirror instability should modulate various 
plasma parameters (e.g., density, temperature, and magnetic field strength) by 
about ± 10%, whereas the resonant-particle energy exceeds the electron temper- 
ature by a factor ~ 10. The amplitude Q of the argument of the exponential 
function in (26) should amount roughly to the product of these two dimension- 
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less factors in the case of a resonant-particle instability, thus Q ~ 1. 


In order to proceed with the analysis of the model, we need to determine 

the group velocity of the electrostatic cyclotron waves v in the frame of the 

S 

drifting flux tube. As far as we are able to determine, observation does not 

shed any light on this point, but, considering the spectrum of particle drift 

velocities involved, v is probably bounded by v =» 0 (waves drifting with the 

S & 

+ — * 

cold plasma in the flux tube) on the one hand and by v = r0_<? (waves drifting 

S — 

with the energetic electrons) on the other. We shall investigate both limits 
to show that the spatial modulation is effective in both cases. 

Model A: Waves drifting with energetic electrons (v » r0_ <p). Defining 

S * 

A A 

the dimensionless variable £ = K/TrO^C^ * K ) » Che model equations for this 
case can be written in dimensionless form. 


ai/ac - - — — + s (27a) 

l + X 

3(Zn X)/3C - (2 Y* T) [— £ — exp (q sin k§) - lj (27b) 

1 + X 

_ — ie — * _ _ * * 

where X = Xt, I =• I/1q, S = tS/Iq anc * k = k t • <) . Numerical solution 

of (27) as a nonlinear initial value problem is implemented. The numerical 
steady state exhibits spatial modulation of the trapped electron flux I, but, 
more importantly, very deep spatial modulations of the precipitation flux 

XI/ (1 + X) is found. Figure 13 shows the solutions for a case with typical 

_ — * 

parameters: S ■ 0.1, Q = 1.0, k * 3.0 and 2 YqT =* 30.0. The scale length of 

A A 

the spatial modulation (furrows) for this case is = (2ir/k) [irfi^ * *0 J • 

Model 3: Waves drifting with cold plasma (v =0). We obtain for this 

S 

case 
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NORMALIZED RATES 



Fig. 13. Electron flux modulation by mirror instability for the case in 
which electrostatic cyclotron waves drift with energetic 
electrons (Model A). 
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31 / 3 ? 


(28a) 


I I 

1 + X 


+ S 


1+1=1 exp(Q sin k ?) 


(28b) 


The algebraic equation (28b) can be solved for X. Because the diffusion rate 
X is restricted to positive or zero values, we have 


X = Max II exp (Q sin k 5) - 1; 0] 


(29) 


whence 


= S - Max [I - exp(- Q sin k?); Oj (30) 

Equation (30) can be solved numerically in the same way as (27). For the same 
parameter values given in Model A, above, the solution of (30) is shown in Fig. 
14. The properties and the scaling of this solution is the same as that of 
Fig. 13 except that, as the result of (29), the precipitation rate at minimum 
is exactly zero. The solution to (30) relaxes to the steady state at a rate 
slower than that of (27); however, this is not of any significance since only 
the steady rate is to be considered in the context of our theory. 

With regard to (30), it is interesting to note that one can seek an 
analytic solution by an expansion of the form 

00 

I = S + I n (Q) + E [C cos nk? + S sin nk?J (31) 

u , n n 

n=l 

by taking advantage of the expansion [e.g. Olver, 1964; p. 37 6 J 
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NORMALIZED RATES 



Fig. 14. Electron flux modulation by mirror instability for the case in 
which electrostatic cyclotron waves drift with cold plasma 
(Model B). 
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exp (-Q sin k£) = I q (Q) + 2 2 (~l) q (Q) cos (2qk5) 

q“i q 

00 

-21 (-l) q I 2 (Q) sin C 2q k? + kC) (32) 

q=0 q 

The results of this exercise confirm the properties of the numerical solution. 
In both models, the scale length of the modulation is 

k q =■ |(2ir/k)[T r 5 3 (♦ • k)J| = 2r/k (33) 

To compare <q with the ion gyroradius scale of a typical 3-keV proton ~ 

km, we note that 8^ T ~ 10”® L® [e.g., Schulz and Lanzerotti, 1974; p. 21; 

A A 

Schulz, 197 4b ] . The angle < • $ can be obtained from (1) by requiring that 

A A 

k ■ i E for reasons documented above. Thus we obtain 

|(r/p Q )(< '• h fl 3 t| ~ 10 " 2 L 6 [L 6 (l + 3 sin 2 ^>) 

3 2 

- 2(4 /$ ) L J sin<? + (<5 /$ )^J cos $ (34) 

r c r c 

For auroral conditions, » (4 r /® c ), so (34) implies 

r fl 3 t (< • 9) ~ 10 -2 P Q L 3 cos * ~ icf 2 L 6 cos <j> km (35) 

Using the assumed values k ■ 3 and L = 7, (33) and (35) imply that the value 
of mapped to the ionosphere is about 6 km. This is roughly in agreement 
with the observed scale modulation of morningside diffusive auroras [Akasofu, 
1 974 J and with patches [Yamamoto and Oguti, 1982J. Further confirmation of 
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spatial modulation at the ~ 10 km scale has recently been obtained [Vondrak et 
al., 1983] by incoherent-scatter radar observations. These were taken 
equatorvard of the diffuse aurora in the evening sector. Although our theory 
prefers modulation in the morning sector because of the observed pattern of 
ion pitch-angle anisotropies obtained by SCATHA, all other elements of the 
theory can be operative at any local time provided that the mirror instability 
can be driven. 

Finally we make a fairly obvious point related to the time dependence of 
(23), which we ignored in the subsequent development but need not be in 
general. In general, the colder convecting plasma can establish the spatial 
modulation of the magnetic field in the manner considered prior to this 
section. The source of the hot electrons, which precipitate and drive the 
aurora, need not be time-independent. Indeed, auroral brightenings associated 
with injection events would imply that the source of hot electrons is defi- 
nitely time-dependent. If so, the mechanisms of spatial modulation can be 
translated into both time and space modulations, as required by (23). The 
natural consequence would be a theory of morningside pulsating patches and 
diffuse auroral striations [e.g. Berkey, 1980] (as distinguished from the 
theory of the rare pulsating arcs [Cornwall and Chiu, 1982J). The oscillatory 
time-dependent structure of (22) , in isolation from the spatial modulation of 
the colder plasma by the mirror instability, has been treated by Davidson 
[1979] who considered a single component of electrons only and, thus, is not 
able to account for the spatial structures and their drifts. According to our 
theoretical mechanism, including time-dependence, pulsating structures (stri- 
ations and patches) would prefer the morningside and, more importantly, would 
not necessarily be associated with micropulsations. Obviously, pulsating 
auroras have local as well as global morphologies [Royrvik and Davis, 1977]. 
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Our observations here are to be taken in the context of a new avenue of 


thought on morningside pulsating auroras, not a global model as yet. However, 
salient observed features such as recurrent patterns [Oguti, 1976 J are easily 
understood in terms of quasi-periodic injections of hot electrons which drift 
over and are caused to precipitate by slowly varying spatial patterns of 
furrows set-up by the warm plasma convection. 
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VII. CONCLUSIONS 


In the above we have constructed a theoretical framework which links 
together diverse observations such as the local time structure of plasma 
characteristics, the global structure of magnetospheric convection, and the 
mechanisms of wave-particle interactions into a coherent picture which offers 
a possible explanation of the morphology of the aurora in the post-midnight 
and morning sectors. This picture complements, indeed demands, the existence 
of upward ionospheric ion acceleration in the evening discrete arc. Thus, in 
a real sense, it embodies the important principle that the aurora arises as 
result of electrodynamic coupling between the magnetosphere and the iono- 
sphere. 

Basing our theory upon observations of auroral and subauroral distribu- 
tion functions by the SCATHA satellite near the equatorial plane, we offer 
kinetic explanations of the morningside auroral morphology. We argue that 
zero frequency mirror instability can be driven by warm ions convected in from 
the plasmasheet and that the instability is stabilized on the eveningside by 
field-aligned ions from the discrete arc. We further argue that these spatial 
modulations of the magnetic field and warm plasma density can modulate the 
precipitation rates of hot electrons according to a morphology selected by the 
action of shear in the convecting warm plasma flow. An essential ingredient 
of our theory is the different actions of the warm and hot plasmas, accounting 
for the observation that morningside auroral structures drift with the warm 
convecting plasma but not with the hot electrons that produce the aurora 
itself. In the present state, our theory represents an assembly of coherently 
related theoretical and observational segments; as such, it is not yet a 
numerical model. A numerical model will be constructed in the future in the 
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context of global arc formation theory integrating the present work with work 
on the evening aurora [e.g. Chiu et al., 1981]. 
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